Effect of high-/^ dielectrics on charge transport in graphene 
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The effect of various dielectrics on charge mobility in single layer graphene is investigated. By 
calculating the remote optical phonon scattering arising from the polar substrates, and combining it 
with their effect on Coulombic impurity scattering, a comprehensive picture of the effect of dielectrics 
on charge transport in graphene emerges. It is found that though high-/^ dielectrics can strongly 
reduce Coulombic scattering by dielectric screening, scattering from surface phonon modes arising 
from them wash out this advantage. By comparing the room-temperature transport properties with 
narrow-bandgap III-V semiconductors, strategies to improve the mobility in single layer graphene 
are outlined. 

PACS numbers: 72.b-g 



Graphene, a 2D gapless semiconductor with a hon- 
eycomb crystal structure of carbon atoms, has gained 
significant attention recently owing to its conical band- 
structure [1], unconventional non-integer quantum Hall 
effect , possible applications emerging from high room- 
temperature (RT) electron mobility (~ 10^ cm^/V-s) |3] 
,and tunable bandgaps in nanostructures carved from it 
in |5] . Appreciable modulation of current by electrostatic 
gating indicate graphene as an excellent material for con- 
ventional and possibly novel electronic devices [6 . High 
mobilities can facilitate coherent ballistic transport over 
large length scales. This can enable novel devices that ex- 
ploit the wavelike features of electrons, resulting in 'elec- 
tron optics' in the nanoscale material [7 . 

For graphene located in close proximity to dielectric 
substrates, the highest RT carrier mobilities experimen- 
tally reported are in the ~ 10^ cm^/V-s range [8j for 2D 
carrier concentrations in the ~ 10^^ /cm^ regime. How- 
ever, suspending graphene by removing the underlying 
substrate and driving off impurities sticking to it resulted 
in much higher mobilities near RT, albeit at low car- 
rier concentrations [9l [TOl [11]. When compared to 2D 
electron gases (2DEGs) in narrow-bandgap III-V semi- 
conductor heterostructures such as In As and InSb [12], 
mobilities and low-field conductivities (charge x mobil- 
ity) in 2D graphene on dielectric substrates reported to 
date are lower. 

Recently, it was predicted that for thin semiconducting 
membranes, impurity scattering can be reduced by sur- 
rounding the membrane with high-Ai: dielectrics [13], and 
a recent experimental observation was made for graphene 
[M]. The reduction of impurity scattering is affected by 
the reduction of Coulomb scattering by the high-A>: di- 
electric. However, several recent works have pointed out 
that high-/^ dielectrics in close proximity with a conduct- 
ing channel in a semiconductor lead to enhanced surface- 
optical (SO) phonon scattering due to remote optical 
phonon coupling between electrons in the channel and 
polar vibrations in the dielectric. This feature presents 
itself not only in graphene [H |T5[ [16] , but in carbon nan- 
otubes [17] , as well as in the workhorse of the electronics 
industry - Silicon Metal- Oxide Field- Effect Transistors 



(MOSFETs) [H]. 

At low electric fields, intra- and intervalley optical 
phonon scattering is negligible due to the insufficient 
kinetic energy of carriers. If a nominal charged im- 
purity concentration of riimp ^ 10^^/cm^ is present at 
a graphene/high-/^ dielectric interface, for carrier con- 
centrations < 10^^ /cm^, acoustic phonon scattering at 
RT is also relatively insignificant owing to the low den- 
sity of states of graphene at the Fermi energy (see 
[IS [201 [m [22]). In this technologically relevant regime, 
the competing effects of impurity and SO-phonon scat- 
tering are responsible for the low-field transport proper- 
ties of graphene. High-/^ gate oxides result in good elec- 
trostatic gate control of the Fermi level in the graphene 
layer, and at the same time reduce Coulombic scatter- 
ing from charged impurities. However, at the same time, 
they give rise to SO phonon scattering. It is impera- 
tive at this stage that this twofold role of dielectrics be 
studied in detail to clarify the competing roles, and con- 
sequently outline routes for improving the mobility and 
mean free path of carriers in graphene. That is the goal 
of this work. 

Surface modes in the phonon spectrum of a solid arise 
due to its finite size. In the case of a polar material, 
the phonon field due to longitudinal surface modes prop- 
agates along the surface and induces a non-vanishing 
decaying electric field outside it [23]. The dispersion 
of these surface optical (SO) phonons and the induced 
electric field in the inversion layer of a semiconductor- 
oxide interface was first calculated by Wang and Mahan 
[24 . Later, Fischetti et. al [18^ solved the dynamical 
dielectric response of a coupled channel/insulator /gate 
system and showed that the channel and gate plasmons 
can alter the dispersion relation of SO-phonon modes at 
semiconductor-oxide interface for a Si-MOSFET struc- 
ture. In this work, we consider the gate as an ideal metal, 
implying that all electric field lines originating from time- 
dependent Coulombic fluctuations at the gate/oxide in- 
terface should terminate on the gate metal [25| produc- 
ing an insignificant effect on the channel (in our case 
graphene) dielectric response. 

Consider the structure as shown in Fig[l| A dielectric 
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FIG. 1: Graphene/Oxide/Metal layer structure, and a 
schematic sketch of the surface-optical (SO) phonon mode 
strength. 



of thickness tox is deposited on a metal gate, and a single 
layer graphene is placed at an equilibrium distance d from 
the oxide, and air {k, = 1) covers the rest of the space. 
Denoting q and r as the two-dimensional wave vector 
and spatial vector in the graphene plane, we write the 
time-dependent SO-phonon field at a point (r, z) as 



(1) 



where the components of SO-modes (j)q^u{z) are solu- 
tions of Maxwell's equation given by [18j 



0, — OO < Z < —tox 

^q,uj{z) = ( 26g,^sinh(7(z + toa:), -tox < z <0 (2) 



dn 



z>0. 



Here, bq^^j and dq^^ are normalization constants. Ap- 
plying electrostatic boundary conditions at the interface, 
we get the secular equation 



+ es{q,uj)tanh{qtox) = 0, 



(3) 



where, Coxi^) and Csiq^oo) are the dynamic dielectric 
functions of the oxide and graphene respectively. For 
bulk dielectrics, the frequency dependent dielectric func- 
tion can be written as 
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Here, e^o is the high frequency (a; ^ oo) permittivity, 
uJj^Q is the n^^ bulk transverse optical phonon frequency 
and fn is the corresponding oscillator strength of the di- 
electric oxide. The values of oscillator strength for differ- 
ent dielectrics are extracted from experimental data as 
outlined in ref. [18 ; the material constants are listed in 
Table H 

In the long- wavelength limit {q 0), the dielectric 
function of graphene in the RPA approximation is [27] 
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FIG. 2: Energy dispersion uj{q) of surface-optical (SO) modes 
for Graphene on Si02 and Zr02. The two SO phonon modes 
labeled ujso,i^^so,2 are nearly dispersionless (except near q ^ 
0), and their magnitudes are close to the values in the absence 
of graphene. The graphene plasmon mode is energetically 
lower than surface modes. 
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Here, = (e^ + l)/2, ujp{q) = VF\/^q^/7rn is 
the plasma frequency of single layer graphene, a = 
e'^ /27r Co hZav^^F is the effective fine structure constant of 
graphene, e is the electron charge, eo is the permittivity 
of vacuum, h = h/27T is the reduced Planck's constant, 
vp ~ 10^ m/s is the Fermi velocity characterizing the 
bandstructure of graphene, n is the 2D carrier density, 
and Kavg is the average dielectric constant of the system. 
For example, for graphene sandwiched between Si02 
{K.si = 3.9) and air {k. = 1), K.avg = (3.9 + l)/2 = 2.45. 
The numerical solution of Eq. [3] has three roots. Two of 
them represent interfacial phonon modes corresponding 
to two bulk TO modes of the dielectric, and the third 
root represents the graphene plasmon mode. Figs [2] (a) 
and (b) show the calculated dispersions of these three 
modes for a graphene layer on Si02 and Zr02. The SO 
phonon modes are found to be nearly dispersionless, ex- 
cept at the lowest wavevectors. The calculation of these 
dispersions show that the SO-phonon modes of the oxide- 
graphene-air system are not too different from the oxide- 
air surface. The effect of the graphene dielectric response 
on the SO modes is negligible, since charges in graphene 
merely screen the electric field lines originating from the 
oxide-graphene interface. From the results in Fig|2j we 
observe that cOpiasmon « ^50,1,^50,2- Hence, for the 
rest of the work we use the calculated dispersionless part 
of the SO phonon modes - the energies of these two modes 
for some of the commonly used dielectric gates are given 
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TABLE I: Surface-optical phonon modes for different dielec- 
tric gates. Parameters have been taken from ref. [18] . 





Si02 


AIN 


AI2O3 


HfOs 


ZrOs 


SiC" 




3.9 


9.14 


12.53 


22.0 


24.0 


9.7 




2.5 


4.8 


3.2 


5.03 


4.0 


6.5 


ojto,i 


55.6 


81.4 


48.18 


12.4 


16.67 




OJTO,2 


138.1 


88.5 


71.41 


48.35 


57.7 






59.98 


83.60 


55.01 


19.42 


25.02 


116^ 


(^SO,2 


146.51 


104.96 


94.29 


52.87 


70.80 


167.58 



^Refererence 15 . 

^Single surface phonon mode measured in reference [26| . 



in Table [D 

With the calculated SO-phonon frequencies above, the 
Hamiltonian of the electron- SO phonon system is given 

by 



int-) 



(6) 



where £{k) = :^hvF\k\ is the kinetic energy of massless 
Dirac fermions close to the (/C, /C') points of the graphene 
Brillouin Zone, and c\{ck) is the electron creation (anni- 
hilation) operator with wavevector k. The second term 
of Eq. (p| is the phonon part of the total Hamiltonian, 



where a^(ap represents the surface phonon creation 
(annihilation) operator for the v^^ mode. The electron- 
phonon interaction part of the Hamiltonian is given by 



the electron wavefunction is assumed to be of the form 
|(z|xe)P = ^{z — (i), is the equilibrium distance of the 
graphene sheet from the oxide surface. In one-phonon 
scattering process, the scattering matrix element is 

(10) 

where, q = |kf — ki| is the difference between the ini- 
tial and final momentum state of electron and 9 is the 
scattering angle. 

The electron- SO phonon scattering rate in the relax- 
ation time approximation [30l [31] is then given by 
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X ^nq^j,5{Skf - Ski - ^^so) + 

(l + n,,,)(5(^fc, (11) 

where n^^j^ = 1/ {ex.-p{hu;q^jy / kT) — 1) is the equilibrium 
phonon occupation number and e2D(g, 0) is Thomas- 
Fermi screening factor of the 2D carriers [32]. Defining 
dimensionless variables u = dq^^y and x = k/q^^y^ where 
/^F is a constant surface-phonon wave- vector, 
the scattering rate can be written in a more compact 
form as 



where the electron-phonon coupling parameter JT^ is 
given by m [H [29] 



'J 7/ 



2Ac 



+ 1 
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(8) 



where, i^^{i^ox) ^^S^ (low) frequency dielectric 

constant of the dielectric, and A is the area of graphene. 
We treat this interaction term perturbatively to calcu- 
late the electron- SO phonon scattering rates. The wave- 
function for the coupled electron-phonon system can be 
written as |^) = |xe)|0)7 where |0) is the phonon vacuum 
state and \xe) is the electron wavefunction of graphene 
close to the Dirac points. The real-space representation 
of the 2D electron wavefunction of graphene near the 
(/C, /C') points is given by the spinor form 



Xe(k,r) = (r|xe) 



ikr 



/2A 



1 



(9) 



where (j) = tdin~ -^{ky/kx), and k, r are the wavevector 
and spatial coordinate in the graphene plane. Since car- 
riers are confined to the graphene plane, the 2:— extent of 
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where, I^{x^u) is a dimensionless integral that can 
be evaluated numerically [33], and ± stands for phonon 
emission (n+^ = 1 + ^q,u) and absorption (n~^ = riq^y) 
respectively. The total scattering rate is obtained by 
summing the emission + absorption processes over all 
SO phonon modes - 



Tphik.T) 



= E 



e[£{k) - fkos 



(13) 



where 6[...] is the Heaviside unit-step function. Fig- 
ure |3ja) shows the SO phonon scattering rates at room 
temperature calculated using the above formalism for 
graphene located on five different polar gate dielectrics. 
For this plot, a 2D carrier concentration of n = 10^^ /cm^ 
and an equilibrium graphene-oxide distance of d = 0.4 
nm is used [15 . Figure |3jb) shows the dependence of the 
SO phonon scattering rate on the graphene-dielectric dis- 
tance for three different values of d. As expected from the 
decay of the SO phonon evanescent mode from the dielec- 
tric surface, the scattering rate reduces with increasing 
d. The kinks in the scattering rates in Fig [3]^ a) indicate 
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FIG. 3: Left - (a): SO phonon absorption + emission scatter- 
ing rates for graphene on various dielectrics. The higher the 
static dielectric constants, the higher the SO phonon scatter- 
ing rate. The onsets of SO phonon emission into two modes 
are indicated by arrows. Right - (b): SO phonon scattering 
rate for graphene on SiO 2 for three different graphene-oxide 
thicknesses. 



the onsets of SO phonon emission processes, as indicated 
by arrows. Zr02 with the highest low- frequency dielec- 
tric constant (^n^ ~ 24) shows the strongest electron-SO 
phonon scattering among the five, and Si02 shows the 
lowest. If SO phonons were the sole scattering mech- 
anism responsible for limiting the low-field mobility of 
graphene, this analysis would indicate that using a low- 
K dielectric would be the most beneficial. However, the 
presence of charged impurities and Coulombic scattering 
due to them changes this simplistic picture. 

The scattering rate due to charged impurities present 
in the graphene-dielectric interface is given by 



,{k) 



m 



WF{a) 



(14) 

where riimp is the sheet density of impurities at the in- 
terface, a = aikp/k) is a dimensionless argument with 
a the effective fine-structure constant defined earlier fol- 
lowing Eq. [5j kp is the Fermi wave vector, and the di- 
mensionless function F{a) is defined in ref. [22]. As is 
clear from Eq. [T4j the impurity scattering rate can be 
strongly suppressed by using high-/^ dielectric materials 
next to graphene. This was pointed out earlier in [13 for 
the general case of semiconductor nanoscale membranes. 
However, using a high-/^ material also enhances surface- 
optical phonon scattering rates, as based on the earlier 
analysis in this work. With increasing dielectric con- 
stant of the graphene environment, impurity scattering is 
damped, but SO phonon scattering increases, indicating 
an optimal choice of the dielectric exists for obtaining the 
highest mobilities for graphene located on substrates. 
Figure [4]^ a) illustrates this fact for graphene on five 



different dielectrics with increasing low-frequency dielec- 
tric constants at room temperature for an impurity den- 
sity of riimp = 5 X 10^^/cm^ and a carrier concentra- 
tion n = 10^^ /cm^. The mobility (/i = a/ne) is calcu- 
lated from the conductivity using the relation cr(T) = 
{e^ /h)2vFkFT. The solid curve is the predicted en- 
hancement of mobility due to the damping of Coulom- 
bic impurity scattering, and the hollow circles indicate 
the expected mobility for the particular dielectrics if SO 
phonon scattering was absent (these are calculated with 
r = Timp)' The filled circles show the degradation of the 
mobility due to SO phonon scattering; for these values 
of mobility, not = \t~^ + ^^mp)'^ is used. It is evi- 
dent that SO phonon scattering nearly washes out the 
improvement of mobility that could have resulted from 
the reduction of impurity scattering, and RT mobilities 
are limited to (Irt ^ 10^ cm^/V-s, as indicated by the 
shaded band. The effect of SO phonon scattering might 
seem minimal for Si02, but it is important to note that 
the SO phonon scattering is independent of the impurity 
concentration; therefore a purer graphene/oxide interface 
will still result in substantial SO phonon scattering and 
mobilities in this regime. In Fig [4]^b) we plot the car- 
rier conductivity in graphene for two different materials, 
Si02 and Zr02. The solid lines correspond to conduc- 
tivity due to impurity scattering alone, while the dashed 
lines incorporate the effect of both impurity and surface- 
optical phonons. The conductivity is linear with carrier 
concentration in agreement with experiments [8 . Due to 
the importance of SO phonon scattering, is also expected 
to be temperature-dependent - with lowering of temper- 
ature, the reduction of SO phonon scattering will result 
in an increase in mobility from the filled circles in Fig 
|4ja) to the hollow circles. 

Among the dielectrics considered here, AIN is the most 
promising for achieving high mobilities due to the higher 
SO phonon frequencies compared to AI2O3, Hf02 and 
Zr02 resulting in lower SO phonon scattering, and a 
higher static dielectric constant than Si02 resulting in 
lower impurity scattering. Along the same lines, epitax- 
ial single-layer graphene grown on SiC substrates are ex- 
pected to be less affected by SO phonon scattering due 
to the high phonon energies of crystalline SiC (see Ta- 
ble [l| also pointed out in [15 ), and at the same time 
have a lower sensitivity to impurity scattering due to the 
higher dielectric constant of SiC compared to Si02 . If one 
compares carrier mobilities in graphene with 2DEGs in 
III-V semiconductor heterostructures, the impurity scat- 
tering component in InAs or InSb based heterostructures 
is strongly damped by modulation-doping, which allows 
the spatial separation of charged impurities from mobile 
carriers. At the same time, narrow-bandgap semicon- 
ductors have intrinsically high static dielectric constants 
{^inAs ^ 12.5 and ninSh ^ 18). Therefore, impurity 
scattering can be effectively damped in such materials, 
and the RT mobilities are typically limited by interface 
roughness scattering (for thin wells), or polar optical 
phonon scattering, which is strong in III-V semiconduc- 
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FIG. 4: (Left - (a): Electron mobility in graphene as a func- 
tion of the gate dielectric constant. High-Ai: dielectrics reduce 
Coulombic impurity scattering, but strong SO phonon scat- 
tering by them reduces the RT mobilities to ^ 10^ cm^/V-s. 
Right - (b): Electron and hole conductivity as a function of 
carrier concentration for graphene on Si02 and Z1O2. 



tors. The fact that graphene is an atomically thin layer 
makes its transport properties especially sensitive to the 
surrounding dielectrics, and by judicious choice of such 
dielectrics, low- field charge mobilities can be improved. 
The ideal dielectrics would be those that possess both 
high static dielectric constants, and high phonon ener- 
gies that are not activated in low- field transport. Among 
the dielectric materials considered, SO phonon scattering 
from them limits the mobilities to lower values than cor- 
responding III-V semiconductors such as In As and InSb. 
Suspending graphene removes the SO phonon scattering 
component, but reduces the electrostatic control by a ca- 
pacitive gate and introduces electro- mechanical effects. 
While this can be a feasible route towards ballistic trans- 
port in graphene, a high-/^ dielectric with high SO phonon 
energies would be the most desirable in the future. 

In conclusion, we have investigated the effect of var- 
ious dielectrics on the electron mobility in single layer 
graphene by considering the effects of impurity and SO 
phonon scattering. By calculating the remote polar opti- 
cal phonon scattering arising from the polar substrates, 
and combining it with their effect on Coulombic impu- 
rity scattering, a comprehensive picture of the effect of 
dielectrics on charge transport in graphene emerges. 
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